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Abstract 

Efimov states are a sequence of sfiallow 3-body bound states tliat arise wfien tlie 2-body scattering 
lengtli is large. Efimov showed that the binding energies of these states can be calculated in terms 
of the scattering length and a 3-body parameter by solving a transcendental equation involving a 
universal function of one variable. We calculate this universal function using effective field theory 
and use it to describe the three-body system of ^He atoms. We also extend Efimov's theory to 
include the effects of deep 2-body bound states, which give widths to the Efimov states. 
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The interactions of nonrelativistic particles (such as atoms) with short-range interactions 
at extremely low energies are determined primarily by their S'-wave scattering length a. 
If \a\ is much larger than the characteristic range / of their interaction, low-energy atoms 
exhibit universal properties that are insensitive to the details of the interaction potential. 
In the 2-body sector, the universal properties are simple and familiar. The differential 
cross section for two identical bosons with relative wave number k <^ l/l and mass m is 
da/dQ = Aa? /{I + k'^o?). If a > 0, there is also a shallow 2-body bound state (the dimer) 
with binding energy B2 = 'h?' Imc? . In the 3-body sector, there are also universal properties 
that were first deduced by Efimov [ij. The most remarkable is the existence of a sequence 
of 3-body bound states with binding energies geometrically spaced in the interval between 

jmc? and Ti' jmt^ . The number of these "Efimov states" is roughly ln(|a|//)/7r if |a| is large 
enough. In the limit \a\ — > cxd, there is an accumulation of infinitely many 3-body bound 
states at threshold (the "Efimov effect"). The knowledge of the Efimov binding energies is 
essential for understanding the energy- dependence of low-energy 3-body observables. For 
example, Efimov states can have dramatic effects on atom-dimer scattering if a > 0, 01 
and on 3-body recombination if a < ^S, lAj • 

A large 2-body scattering length can be obtained by fine-tuning a parameter in the 
interatomic potential to bring a real or virtual 2-body bound state close to the 2-atom 
threshold. The fine tuning can be provided accidentally by nature. An example is the '^He 
atom, whose scattering length a = 104 A is much larger than the effective range / ~ 7 
A. Another example is the 2-nucleon system in the ^5*1 channel, for which the deuteron 
is the shallow bound state. This system provided the original motivation for Efimov's 
investigations 0]. In the case of atoms, the fine tuning can also be obtained by tuning 
an external electric field C| or by tuning an external magnetic field to the neighbourhood 
of a Feshbach resonance 7[. Such resonances have, e.g., been observed for ^^Na and ^^Rb 
atoms 0] and are used to tune the interactions in Bose-Einstein condensates. An important 
difference from He is that the interatomic potentials for Na and Rb support many deep 
2-body bound states. 

Efimov derived some powerful constraints on low-energy 3-body observables for systems 
with large scattering length They follow from the approximate scale- invariance at length 
scales R in the region I <^ R <^ \a\ and the conservation of probability. He introduced polar 
variables H and ^ in the plane whose axes are 1/a and the energy variable sgD.{E)\mE\^^'^ /h, 
and showed that low-energy 3-body observables are determined by a few universal functions 
of the angle ^. In particular, the binding energies of the Efimov states are solutions to 
a transcendental equation involving a single universal function of ^ [^. In this paper, 
we calculate this universal function for the case of 3 identical bosons and apply Efimov's 
equation to the ^He trimer. We also extend Efimov theory to atoms with deep 2-body bound 
states. 

The existence of Efimov states can easily be understood in terms of the equation for 
the radial wave function f{R) in the adiabatic hyperspherical representation of the 3-body 
problem The hyperspherical radius for 3 identical atoms with coordinates ri, r2, 

and r3 is R^ = {rf2 + '"is + '"is)/^) where Vij = \ri — rj\. If |a| ^ /, the radial equation for 3 
atoms with total angular momentum zero reduces in the region I <^ R <^ \a\ to 
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f{R) = EfiR), (1) 



where Sq ~ 1.00624. This looks like the Schrodinger equation for a particle in a one- 
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dimensional, scale-invariant l/E? potential. If we impose a boundary condition on f{R) at 
short-distances of order /, the radial equation ((T)) has solutions at discrete negative values 
of the eigenvalue E = —B3, with ranging from order fi^/mP to order /ma?. The 
corresponding eigenstates are called Efimov states. As \a\ — > cx), their spectrum approaches 
the simple law b'^'' ~ hlh'^h^ / ma? . 

Efimov's constraints can be derived by constructing a solution to Eq. (^J) that is valid 
in the region / -C i? <C |a|. In the case of a bound state with energy E = — -B3, the radial 
variable is = mB^/?? + and the angular variable ^ is 



^ = — arctan(aYmi?3//i) — 7r^^(— a) , (2) 

where 9{x) is the unit step function. Since we are interested in low energies ~ /ma'^, 
the eriergy eigenvalue in (0) can be neglected. The most general solution therefore has the 
form 

/ (R) = VhR [Ae^'« ^""^"^^ + 5e-^'° ^^hr)]^ ^ 

which is the sum of outgoing and incoming hyperspherical waves. The dimensionless coeffi- 
cients A and B can depend on ^. At shorter distances R ~ / and longer distances i? ~ |a|, 
the wavefunction becomes very complicated. Fortunately, we can avoid solving for it by 
using simple considerations based on unitarity [ij. 

We first consider the short-distance region. Efimov assumed implicitly that there are no 
deep 2-body bound states with binding energies B ^ jma? . Thus the 2-body potential 
supports no bound states at all if a < and only the dimer with binding energy B2 = T? /ma^ 
if a > 0. We will address the complication of deep 2-body bound states later. The probability 
in the incoming wave must then be totally refiected by the potential at short distances, so 
we can set B = Ae^^ . The phase 9 can be specified by giving the logarithmic derivative 
Rof'{Ro)/f{Ro) at any point I Rq <ti \a\. The resulting expression for 9 has a simple 
dependence on H: 

9/2 = soHH/cA,). (4) 

The denominator cA* is a complicated function of Rq and -Ro/'(-Ro)//(-Ro)- It differs by an 
unknown constant c from the 3-body parameter A^, introduced in Ref . . 

We next consider large distances ~ |a|. In general, an outgoing hyperspherical wave 
incident on the ~ \a\ region can either be reflected or else transmitted to -R — ^ cxo 
as a 3-atom or atom-dimer scattering state. The reflection and transmission amplitudes 
are described by a dimensionless unitary 3x3 matrix that is a function of only. For 
— TT < ,^ < — 7r/4, scattering states with R —>■ 00 are kinematically not allowed. The 
probability is therefore totally reflected, so we must have B = Ae*^^^\ where the phase 
A depends on the angle ^. Compatibility with the constraint from short distances requires 
^ = A(^) mod 271. Using Eq. ^ for 9 and inserting the expression for H, we obtain Eflmov's 
equation P] 

53 + ^ = ^ e^'^"/'^" exp [A (0 /so] , (5) 
ma^ m 

where ^ is given by Q and the constant c was absorbed into A(^). Note that we measure 
i?3 from the 3-atom threshold and A^, is only deflned up to factors of exp[7r/so]. Once 
the universal function A(,^) is known, the Eflmov binding energies B^ can be calculated 
by solving Eq. (jSJ for different integers n. This equation has an exact discrete scaling 
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FIG. 1: The energy variable — (mi?3/?i^A^)-'^/^ for three shahow Efimov states as a function of 
sgn(a)(A.|a|)~V2. 

symmetry: if there is an Efimov state with binding energy for the parameters a and A^., 
then there is also an Efimov state with binding energy X^B^ for the parameters X~^a and 
A^, if A = exp[n'7r/so] with n' an integer. 

The universal function A(^) could be determined by solving the 3-body equation for the 
Efimov binding energies in various potentials whose scattering lengths are so large that effec- 
tive range corrections are negligible. It can be calculated more easily by using the effective 
field theory (EFT) of Ref. 0] in which the effective range can be set to zero. In Ref. 0, 
the dependence of the binding energy on a and A,, was calculated for the shallowest Efimov 
state and a > 0. In order to extract the universal function A(^), we have calculated the 
binding energies of the three lowest Efimov states for both signs of a. In Fig. ^ we plot 

function of sgn(a)(A*|a|) for these three branches of Efimov states. 
The binding energies for deeper Efimov states and for shallower states near (A^lal)"^/^ = 
can be obtained from the discrete scaling symmetry. A given 2-body potential is character- 
ized by values of a and A,, and corresponds to a vertical line in Fig. ^ such as the dashed 
line shown. The intersections of this line with the binding energy curves correspond to the 
infinitely many Efimov states. Those states with B^ > h'^ /mP are unphysical. For B3 — »• cxo, 
the angle ^ goes to — 7r/2. The ratio of the binding energies of successive Efimov states 
therefore approaches exp[27r/so] ~ 515. However, for the shallowest Efimov states, this ratio 
exhibits significant deviations from the asymptotic value. If a > 0, there is an Efimov state 
at the atom-dimer threshold B^ = h'^/ma? when soln(aA^,) = 1.444 mod vr. The sequence 
of binding energies i?3 in units of Ti^ /ma? is 1, 6.8, 1.4 x 10^, ... . Consequently, the ratio 
of i?3 for the two shallowest Efimov states can range from 6.8 to 210. If a < 0, there is an 
Efimov state at the 3-atom threshold B^ = Q when sq ln(aA*) = 1.378 mod tt. The sequence 
of binding energies -B3 in units of li? /ma? is 0, 1.1 x 10'^, 6.0 x 10^, ... . Thus, the ratio of -B3 
for the two shallowest Efimov states can range from 00 to 550. 
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Using Eq. (0) with n = 0, we have extracted the universal function A(^) from the data 
for the middle branch in Fig. ^ The extracted values of A (^) are given in Table HI An 
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-0.933 
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TABLE I: The values of the universal function A(^). 



analytic expression for A(,^) is not known, but we can obtain parametrizations in various 
regions of ^ by fitting the data: 
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: A = 3.10x2 - 9.63a; - 2.18, 



8. -A 

f]:A = 



= 1.17^ + 1-97^ + 2.12?/ - 8.22, 
0.25^2 + 0.282-9.11, 



(6) 
(7) 
(8) 



where x = (— 7r/4 — ^)^/^, y = ti/2 + ^, and z = {jr + ^y exp[— l/(7r + ,^)2]. These parametriza- 
tions deviate from the numerical results in Table U by less than 0.013. The discontinuity 
at ^ = — ^ and ^ = — is less than 0.016. This accuracy is sufficient for most practical 
calculations using Eq. (0). 

Universality can be exploited to greatly reduce the calculational effort required to predict 
3-body observables for atoms with large a. The observables can be calculated once and for 
all as functions of a and A,, either by using EFT or by solving the Schrodinger or Faddeev 
equations with various methods. The binding energies obtained by solving Efimov's equation 
(0) are shown in Fig. |21 Simple expressions can be given for other observables, such as the 
S-wave atom-dimer scattering length: 

ai2 = a (1.46 - 2.15 tan[so ln(aA,) + 0.09]) , (9) 

as well as the phase shifts and the rate constant for 3-body recombination at threshold • 
Given a and a measured or calculated value of B3 for one Efimov state as input, one can 
read off A^, from Fig. El Predictions for other 3-body observables, such as the atom-dimer 
scattering length in (jU)), are then immediate. 

One of the most promising systems for observing Efimov states is ^He atoms. The ^He 
trimer has been observed 3| , but no quantitative experimental information about its bind- 
ing energy is available to date. The binding energy has been calculated accurately for various 
model potentials. They indicate that there are two trimers, a ground state with binding 
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FIG. 2: The Efimov binding energies Bsma^/h^ as a function of |a|A* for a > (solid lines) and 
a < (dashed lines). Vertical dashed line gives |a|A,, for LM2M2/TTY potentials. Horizontal 
dotted line is the atom-dimer threshold (a > 0). 



energy B^^ and an excited state with binding energy 



The most accurate calculations 



have been obtained by solving the Faddeev equations in the hyperspherical representation 
Till, in configuration space and with hard-core boundary conditions 3|- These meth- 



ods all give consistent results. The results of Ref. [15[ for B2, Bf\ B^'' 
scattering length for four different model potentials are given in Table |n] The discrep- 



i?3 , and the atom-dimer 
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B2 


^(0) 


^(1) 


ai2 


as A* 


^(0) 

^3 


«12 


HFDHE2 


0.830 


116.7 


1.67 




1.258 


118.5 


87.9 


HFD-B 


1.685 


132.5 


2.74 


135(5) 


0.922 


137.5 


120.2 


LM2M2 


1.303 


125.9 


2.28 


131(5) 


1.033 


130.3 


113.1 


TTY 


1.310 


125.8 


2.28 


131(5) 


1.025 


129.1 


114.5 



TABLE II: The values of B2, B^^\ -83^^ and 012 for four model potentials from Ref. the value 
of ttB^* determined from .83^^ and the predictions for B^^ from Eq. © and 012 from Eq. ©. All 
energies (lengths) are given in mK (A), {h^/m = 12.12 KA^ for ^He). 



ancies between B2 and the large-a prediction h^/ma? are about 6-8% 3]- They can be 
attributed to effective range corrections and provide estimates of the error associated with 
the large-a approximation. 

We proceed to illustrate the power of universality by applying it to the ^He trimer. For 
the scattering length a, we take the value a_B = h,j\pmB^ obtained from the calculated 
dimer binding energy. We determine A* by demanding that B^^ satisfy © with ra = 1. 
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Solving Eq. (jSJ with n = 2, we obtain the predictions for S3 in the second-to-last column 
of Table ITTl The predictions are only 1-4% higher than the calculated values, which is within 
the expected error for the large-a approximation. This demonstrates that the ground state 
of the ^He trimer can be described by Efimov's equation (jSJ. If we use the calculated values 
of a as input instead of B2, the predicted values of B^^ are larger than the calculated values 



We can use the value of a^A^ determined from the excited state of the trimer to predict 
the atom-dimer scattering length ai2 and compare with the calculated values in the fourth 
column of Table ITTl The predictions for the four model potentials are given in the last column 
of Table |n] They are smaller than the calculated values by about 13%. If the calculated 
value of a is used as input they are smaller by about 28%. It should be possible to account 
for these differences quantitatively by taking into account higher order corrections fll']. 




Efimov implicitly assumed in his derivation of Eq. (jSJ that there were no deep 2-body 
bound states [H. If such states are present, the Efimov states become resonances that can 
decay into a deep 2-body bound state and a recoiling atom. Thus their energies are given 
by complex numbers E = — -B3 — If a potential supports many 2-body bound states, 

the direct calculation of the widths by solving the Schrodinger equation is very difficult 
[16^ . However, one can show that the cumulative effect of all deep 2-body bound states 
on low-energy 3-body observables can be taken into account by including one additional 
low-energy parameter rj. Three-body recombination into a deep 2-body bound state with 
binding energy can only take place if i? ~ /. It is obvious that the atoms 

that form the bound state must approach to within a distance of order /, since the size 
of the bound state is of order /. However, the third atom must also approach the pair to 



within a distance of order /, because it must recoil with momentum y4mi?/3 ~ h/l, and 
the necessary momentum kick can be delivered only if i? ~ Z. The atoms can approach such 
short distances only by following the lowest continuum adiabatic hyperspherical potential, 
because it is attractive in the region / -C i? -C \a\, while all other potentials are repulsive. 
Thus all pathways to final states including a deep 2-body bound state must fiow through 
this lowest continuum adiabatic potential. Note the lowest continuum state can either be a 
3-atom state or an atom-dimer state. See, e.g.. Fig. 5 in Ref. ^3] for an explicit calculation 
of these potentials. The cumulative effects of 3-body recombination into deep 2-body bound 
states can therefore be described by the refiection probability e"'' for hyperspherical waves 
entering the region i? ~ / of this potential. Up to corrections suppressed by l/\a\, the 
low-energy 3-body observables are all determined by a. A*, and 77. 

We proceed to generalize Efimov's equation (0) to the case in which there are deep 2- 
body bound states. Again, we combine the analytic solution (jSJ to the radial equation in the 
region I <^ R <^ \a\ with simple probability considerations at i? ~ \a\ and R I. Since the 
existence of deep 2-body bound states plays no role in the unitarity constraint at i? ~ |a|, we 
still have B = Ae*^*^^-* with the same universal function A(^) given in (jHl)-®- In the unitarity 
constraint at i? ~ Z, we need to take into account that only a fraction of the probability 
that fiows to short distances is refiected back to long distances through the lowest adiabatic 
hyperspherical potential associated with the 3-atom or atom-dimer continuum. Denoting 
the refiection probability by e~^, the coefficient -B in (jS)) can be written B = Ae'^+'^Z^. The 
compatibility condition 6 — iri/2 = A{C,) mod 2it then becomes 



by 11-21%. 





B3 + -Eg + 



* 27rn/so {A({)+i»7/2)/so 



(10) 
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where ^ is defined by 1^ with i?3 + zr3/2. To solve this equation, we need the analytic 

continuation of A(^) to complex values of ^. The parametrization ©-(IHl) for A(^) should be 
accurate for complex ^ with sufficiently small imaginary parts, except perhaps near ^ = — vr 
where it has an essential singularity. If and for one Efimov state are known, they 
can be used to determine A,, and rj. The remaining spectrum of Efimov states and their 
widths can then be calculated by solving Eq. (jiup . For infinitesimal t], the widths approach 
Fa (?7/so)(-B3 + /ma?') . The widths of the deeper Efimov states increase geometrically 
just like the binding energies, as has been observed in numerical calculations [l6| . If 77 is so 
large that i?3 ~ F3, the Efimov states cease to exist in any meaningful sense. 

We have calculated the universal function A(^) that appears in Efimov's equation 
for the binding energies B3 of Efimov states. This equation can be used as an operational 
definition of the 3-body parameter A^, introduced in Ref. 0]. If B3 for one Efimov state is 
known, it can be used to determine A*, and then universality predicts all low-energy 3-body 
observables as functions of a and A*. In Eq. (Iiup . we have generalized Efimov's equation 
to permit deep 2-body bound states. The generalization involves an additional inelasticity 
parameter t], but the spectrum is determined by the same universal function A(^). 
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